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O^l ' Abstract. Characteristic examples of continuous symmetries in hydrodynamic plasma 

£H , theory (partial differential equations) and in kinetic Vlasov-Maxwcll models (integro- 

& ' differential equations) are considered. Possible symmetry extensions conditional and ex- 

tended symmetries are discussed. Physical reasons for these symmetry extensions are 
r**"* ' clarified. 
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1 Introduction 



Symmetry considerations essentially help us to solve nonlinear problems of plasma physics. 
Among the recent publications, a survey of methods which allow us to find symmetries of integro- 
differential equations of kinetic plasma theory pQ and a review of applications of symmetry 
| methods to the hydrodynamic plasma models based on the partial differential equations [2] can 

be cited. 

In the present paper, some characteristic examples of the symmetries are presented for var- 
ious plasma theory models: electron magnetohydrodynamics, collisionless electron plasma os- 
qq ' dilations, drift waves in plasma and multi-component collisionless plasma containing particles 

■ with equal charge to mass ratios. 

In Section [21 Lie point symmetries are obtained for the partial differential equations of the 
electron magnetohydrodynamics by means of the standard Maple 11 program. A simple invariant 
5_i ■ solution is presented, which is in fact a background solution of a perturbation theory formalism. 



An additional condition is underlined which can lead to the symmetry extension. 

In Section El continuous symmetries of the integro-differential kinetic equations of the colli- 
sionless electron plasma model are discussed. The symmetries were found in [3] by an indirect 
algorithm which allows us to obtain symmetries of the kinetic equations from the symmetries of 
an infinite set of partial differential equations for the moments of distribution functions. 

In Section HJ some conditional symmetries are presented for the Hasegawa-Mima hydrody- 
namic model describing drift waves in a plasma. Due to these symmetries, invariant solutions 
are possible describing the asymptotic structure of the nonlinear waves. 

In Section [51 additional symmetries for integro-differential equations of the kinetic theory 
of collisionless plasma containing particles with equal charge to mass ratio are considered. For 
example, alpha particles and deuterium ions participating in a thermonuclear reaction D + + 
T + — > He ++ + n + 17.6 MeV have close charge to mass ratios. By simple vector considerations 



*This paper is a contribution to the Proceedings of the Seventh International Conference "Symmetry in 
Nonlinear Mathematical Physics" (June 24-30, 2007, Kyiv, Ukraine). The full collection is available at 
|http: / /www.emis.de/journals / SIGMA /symmetry2007.html | 
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so called extended symmetry transformations are obtained which allow us to reduce the number 
of equations. 

Conclusions are made in Section [6J 



2 Electron magnetohydrodynamics 

Let us consider the equations of the electron magnetohydrodynamics (EMHD) [3]: 

— = V x (v x * = B-AB, v = -VxB, V • B = 0, 

at y ' 

where ^ is the generalized vorticity, B is magnetic field strength and v is the hydrodynamic 
velocity of the electron plasma component. Nine functions, if, v, B, of four independent variab- 
les t, x, y, z are present in this system of partial differential equations. Nevertheless, its 
symmetry has been obtained by the standard Maple 11 package program. 

The infinitesimal operators of Lie point symmetries of the model are as follows: 

x - 3 x - 9 x - 3 x - 3 
at ox oy oz 

d 3 d d 

d d „ d T d 

Xq = r x — — h v x — — h B x h * x -— . 

b dr Ov dB 

The Maple 11 program was tested previously on Korteweg-de Vries, nonlinear Schrddinger, 
Hasegawa-Mima and other nonlinear plasma theory models and reproduced well known sym- 
metries. So we can expect that operators X\— Xq form a full basis of EMHD Lie point symmetry 
algebra. 

A little doubt remains, however: maybe some additional symmetries could be obtained by 
the direct 'manual' use of the standard Lie algorithm. 

Due to time and space homogeneity (symmetries X\ to X4) the simplest exact solution exists 
corresponding to the constant external magnetic field which can be directed along the Oz axis 
without the loss of generality (because of the rotation symmetry Xq). The field amplitude can 
be set as Bo = 1 due to the similarity transform generated by X^\ 

B = e z , v = 0. (1) 

This solution can be chosen as a background for the perturbation theory. After the substi- 
tution B^e^ + B, VP — > e z + $ the first EMHD equation becomes 

9* dv 

_ = _ + Vx(vx *) (2) 

but other equations remain unchanged. In this way the system of equations describing helicon 
waves in a plasma is usually obtained. 

We can expect that EMHD model is much more symmetric. For example, nonlinear term 
in ([2]) vanishes under an additional condition: 

vx$ = VF, (3) 

where F is arbitrary scalar function. In this way we obtain the set of linear partial differential 
equations and the nonlinear algebraic restriction This property can be a source of rich 

additional symmetries. 
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3 Collision less electron plasma 

Let us consider electron collisionless plasma. In this case the Vlasov-Maxwell integro-differential 
system of equations holds: 

df_ v df__ df_ _ 
dt dx dv ' 



oc 



dE r°° dE 

^ =1 -j_J dv > -dt=L vfdv ^ (5) 

where f(t,x,v) is the distribution function of the electron component of the plasma, E(t,x) is 
the electric field strength. 

We cannot find symmetries of this integro-differential model by the traditional Lie method. 
This can be done, however, by the indirect algorithm [3]. 

First, let us introduce the moments: 

/■oo 

M k (t,x)= / v k fdv, k = 0,1,2,... (6) 



and obtain an infinite set of partial differential equations for them, which is as follows: 

dM k dM k+l dE „, dE „, 

dt ox ox dt 

Then, we consider the reduced system involving only first N + 1 equations of this set, and 
we can obtain by the standard Lie procedure the symmetries of this system: 

d d „ d „ d , . , d 



k=l 



d d\ d 



X * = C ° S W ( d~x + dE ) ~ Sm W g kM ^8M k ' (7) 

d 



(d d \ jv+i 
~dx~ + ~dE) +cosW E fcMfc 
' k=l 



dM k 



and an additional symmetry 



d ( dF(t)\ d 



X FGN = F(t)——+[G(t)-x 



dM N V dt J dM N+1 

F(t) and G(t) being arbitrary functions of t. 

The standard Maple 11 program can help us to check these results for not very large TV 
(depending on computer capacity). 

Only two highest moments M/v and M^ + \ are present in Xfqn- So only moments with 
numbers greater than any previously chosen are involved in the transform generated by Xfgn- 
In the limit N — > oo the symmetry group is generated by X\— X§. 

In general, restoring kinetic symmetries from the symmetries of the corresponding set of the 
moments equations is a hard problem in pure mathematics. Nevertheless, for the relatively 
simple operators X\-X§ it was shown in [3] that in the case considered they correspond to the 
following infinitesimal operators of the kinetic theory: 
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/ d d \ d ( d d \ d 

X 4 = cos (t) + -J - sin (*)_, X 5 = sin + cos (*)-. (9) 

Strictly speaking, from the definition © it can be readily seen that the transformations gen- 
erated by the operators (jHJ), © in the space of variables t, x, v and / lead to the transformations 
generated by the infinitesimals (|7|) in the space of variables t, x and M^, k = 0,1,2, ... . 

The symmetries of the same electron collisionless plasma kinetic model ([4]), © were found 
in [5] by the direct algorithm, i.e. without using of an infinite set of moments equations. The 
same generators ©, © were obtained. 



4 Hasegawa— Mima model 

Let us consider an inhomogeneous plasma slab in the external homogeneous magnetic field. 
Electrons, unlike ions, are magnetized, smoothing an electrostatic potential $ along the magnetic 
field lines. In this case, Hasegawa-Mima model equations hold [6]: 

_ + J($,*) = _, * = d>-A ± <I>, (10) 



where ^ is the generalized vorticity, 

j(FG) = — — -— — A = — + — 

dx dy dx dy ' dx 2 dy 



The simultaneous presence of the $ and in the second equation of (|1U|) leads to the 

symmetry reduction and, as a consequence, to the absence of self-similar solutions [7]. 

Nevertheless, self-similar solutions can exist as a consequence of conditional symmetries [8]. 
The notion of conditional symmetry is discussed and large bibliography is presented in a recent 
paper [9]. 

For example, if an additional condition 
<9 2 $ 

(11) 



dx 2 



is fulfilled, the term Aj_$ is reduced to which leads to the symmetry extension. In this 
way we obtain more symmetric, but still nonlinear, set of equations, which admits the following 
similarity transform: 

9 d T d 

Xi = x— + <£- h * , 

dx 9$ 

which allows us to try solutions in the form: 
* = xF(t,y), V = xG(t,y). 

G = — 1 is a trivial solution which is not interesting for physics. So only the simple nonlinear 
exact solution with G / -1 important for physics as a nonlinear wave asymptotic, will be 
considered below. More complicated solutions can be presented and discussed in future papers. 

For the functions F(t, y) and G(t, y) we obtain the equations 

+ | d(F/(G + l)) G = F _^£ (12) 

dt dy dy 2 

Equations (|12p have a particular solution [8]: 



F = a(l + /3cos((w + 8uS)t + qy)), 



(13) 
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so that 

$ = ax(l + (3 cos((uj + 5u))t + qy)), (14) 

where a is arbitrary constant amplitude, constant factor j3 determines the relative weight of the 
zonal flow and the monochromatic wave, the frequency is uj = q/(l + q 2 ) and the frequency shift 
is 5lo = -aq 3 /(l + q 2 ). 

Another class of conditionally symmetric invariant solutions became possible under the fol- 
lowing condition: 

d 2 $ _ 
dy 2 

In this case the following similarity transform appears as a conditional symmetry: 

X -t— — 

2 dt ^ dy 

and another class of invariant solutions is possible 

$ = F(x)| + G(x), ^ = (F(x) - F"(x))| + G(x) - G"(x). 

Pure self similar solutions G{x) = are absent in this case, but there exist solutions corre- 
sponding to the arbitrary shear flow F(x) = 0, G(x) being an arbitrary function of x. 

5 Multi-component collisionless plasma 

Let us consider TV-component collisionless plasma. In this case the Vlasov-Maxwell integro- 
differential system of equations holds: 

a/Q +v^ + ^fE + l[vxBA^ = 0, (15) 



dt dr m a \ c J <9v 

VxE + ^ = 0, V-E = 4vrp, VxB = ^ + ^j, V • B = 0, (16) 

where r, v) is the distribution function of the a—th component of the plasma, a = 1, . . . , N. 
Charge and mass of particles of the ath component are denoted by e a and m a , respectively. 
Charge and current densities have the form 



N poo N poo 

p = ^2e a / f a dv, j = ^ e Q / vf a dv. (17) 

a=l J -°° a=l J -°° 

Let a = 1 and a = 2 correspond to plasma components with equal charge to mass ratio of 
particles, for example, alpha particles and deuterium ion components, so we assume (and it is 
a good approximation: 

ei _ £2_ 
mi ni2 

According to (1170 . these components enter the Maxwell equations only as the sum 

ei/i + ea/a, (18) 
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which is in fact the distribution function of the charge density of the components a = 1 and 
a = 2. In addition, according to (|15|) . the functions f\, f 2 and e±fi + 62/2 satisfy the same 
Vlasov equation. So the transformation [8] 

fi = ft ~ e 2 F(h, f 2 ), ti = f 2 + eiF(f u f 2 ), (19) 

where F(fi,f 2 ) is an arbitrary function of its arguments, leaves the Vlasov-Maxwell equations 
invariant, at least if we perform in (|17p the summation first and the integration later. 

The transformation (|19p does not preserve, in general, the positiveness of the distribution 
functions. Moreover, they can lead to the divergence of the moments, which can be shown, for 
example, by choosing F = const. 

Nevertheless, we can choose F(f±,f 2 ) = fi/e 2 to obtain the simplifying transformation 

f[ = 0, e 2 f 2 = e l f l + e 2 h. (20) 

The simplifying transform (|20j) can also be obtained from the symmetry deduced from the 
symmetry of moments' equations, 

fi = h exp(a), f 2 = f 2 + — (1 - exp(a))/i, 

e 2 

a is arbitrary constant, in the limit a — > — 00. 

According to (|20p . we can choose the invariant (|18p as a new distribution function and solve 
the system (fl~5j) . ([IB]) for 

ei/i + e 2 / 2 , /3, In 
and the fields E and B, so the number of equations is reduced by one. 



6 Conclusions 

Lie point symmetries of electron magnetohydrodynamics are obtained (Section [2]), which include 
time and space homogeneity, self similarity and rotations. As a consequence of the symmetries, 
simple invariant solution (Q]) exists which is in fact the background solution which allows to 
explore nonlinear waves evolution by a perturbation formalism. Rich conditional symmetries 
are possible in this model due to additional conditions like ([3]). 

Lie point symmetries of the integro differential Vlasov-Maxwell model for collisionless elec- 
tron plasma obtained in [3] by the moments' method, i.e. from the symmetry of an infinite set 
of partial differential equations for the moments of the distribution functions, are presented in 
Section [3J 

It is shown in Section [4] that hydrodynamic Hasegawa-Mima model for drift waves in magne- 
tized plasma, which, in general, is not self similar, has exact asymptotic self similar solutions (|14p . 
It is possible due to the conditional symmetry of the model as a consequence of an additional 
condition (|lip . The physical reason for this symmetry extension is that nonlinear waves have 
a simpler shape near their critical points - nodes, crests etc. 

In Section [5j by simple vector considerations, it is shown that multi-component collisionless 
plasma containing components with equal charge to mass ratio of particles is invariant under 
the transformation (I19p . if we perform the summation first and the integration later in the 
expressions for charge and current densities (|17p . This transformation, called the extended 
symmetry, does not preserve, in general, the global conditions of non negativity of distribution 
functions and existence of their moments. So, in contrast with the conditional symmetries, this 
symmetry extension is a consequence not of imposing additional conditions, but of neglecting 
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some global constraints. Nevertheless, among the extended symmetries (|19j) we can find the 
simplifying transform (|2U|) . which allows us to reduce the number of equations by 1. 

The physical reason for the extended symmetry (|19p is that particles with equal charge to 
mass ratios under the same initial conditions (coordinate and velocity) in a given electric and 
magnetic fields move along the same trajectories. The second physically important reason is 
that we consider very fast evolution of a plasma, faster than any collisions. 
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